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ABSTRACT 

A simple formalisni to describe nonthermal electron acceleration, evolution, and radiation in super- 
nova remnants (SNRs) is presented. The electron continuity equation is analytically solved assuming 
that the nonthermal electron injection power is proportional to the rate at which the kinetic energy of 
matter swept up in an adiabatically expanding SNR shell. We apply this model to Fermi and HESS 
data from the SNR RX J1713. 7-3946, and find that a one- zone leptonic model with Compton-scattered 
cosmic microwave background (CMB) and interstellar infrared photons has difficulty providing a good 
fit to its spectral energy distribution, provided the source is at a distance ~ 1 kpc from the Earth. 
However, the inclusion of multiple zones, as hinted at by recent Chandra observations, does provide 
a good fit, but requires a second zone of compact knots with magnetic fields B ^ 16 /iG, comparable 
to shock-compressed fields found in the bulk of the remnant. 

Subject headings: supernova remnants — acceleration of particles — radiation mechanisms: nonther- 
mal — Shockwaves — gamma rays: theory 



1. INTRODUCTION 

Acceleration of particles at SNR shocks is considered 
the leading mechanism for the production of cosmic-ray 
protons and ions from energies of ~ 1 GeV/nucleon 
up to the knee of the cosmic-ray spectruni (e.g. , 
Lagage fc Cesarskvl [l983t IBlandford fc Eichleri [l987l: 



Berezinskii et al.l I199O0 . As the expanding SNR shell 



sweeps up matter from the surrounding circumstellar 
medium (CSM), a pair of shocks is formed, with the lead- 
ing forward shock sweeping up and accelerating CSM ma- 
terial, and a reverse shock braking the metal-rich SNR 
ejecta. A few select particles gain energy as they ran- 
domly diffuse back and forth across each of the two shock 
fronts, while convecting downstream into the shocked 
fluid. In the test particle theory of first-order Fermi ac- 
celeration, this leads to a nonthermal power-law particle 
distribution in momentum with number spectral near —2 
for a compression ratio near 4. The spectrum that re- 
sults after folding in the effects of diffusive escape from 
the SNR into intergalactic space, and from the disk of 
the Galaxy into the halo, is in reasonable agreement with 
the ineasured co s mic-ray spec trum fe.g.. Jones fc Ellison! 
[I99l iKnUiggl iHiilasI [20051) . 

Charged cosmic rays are deflected by Galactic mag- 
netic fields during transport, so their direction does not 
point back to the original production site. This is the 
principal reason that the sources of Galactic cosmic rays 
remain elusive a century after their discov ery. Given 
the ex treme difficulty in detecting neutrinos ()Yuan et al.l 
l2011bD . electromagnetic signatures of cosmic rays offer 
at present the best opportunity to identify the sources of 
cosmic rays and settle the question of their origin (e.g., 
Ginzburg fc Svrov atskiilll964l:lGaiss"CTlll990l : lD7'urv et alJ 
~994 : Reynolds 200l). 



With speeds reaching IC* km s~^ or more, SNR shocks 
are likely candidates for accelerating electrons and nu- 
cleons to high energies. The polarization and spec- 
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tral properties of the smooth broadband nonthermal ra- 
dio through X-ray emission is almost certainly electron 
synchrotron radiation, though the optical, UV and X- 
ray spectra can additionally reveal s trong line signa- 
tures from shock- hea ted shell material (|Slane et al.ll2002l : 
iBadenes et al.l[2006( ). Several radiative mechanisms can 
be responsible for 7-ray emission. Energetic leptons emit 
7 rays through Compton scattering of the ambient ra- 
diation fields, principally the cosmic microwave back- 
ground (CMB) and ambient stellar and IR fields, but 
also by scattering photons of the the synchrotron field. 
Electrons make nonthermal bremsstrahlung 7 rays when 
colliding with target gas and dust particles in the cir- 
cumstellar medium (CSM) or in the shocked shell ma- 
terial. Nuclear collisions of hadrons make pion-decay 
7 rays when cosmic-ray protons and ion interact with 
that same matter. Spectral and morphological differ- 
ences are expected between a leptonic and hadronic ori- 
gin. Most well known is the prediction of the 7r° de- 
cay feature peaking at 70 MeV in a photon spectrum 
(jGinzburg fc Svrovatskii|[T964l : lHavakawalll969() that re- 
sults from hadronic processes. For an electron injection 
spectrum softer than number index g = 2 (where the in- 
jected electrons are injected with spectrum (3(7) oc 7""^), 
as expected in the test-particle limit for strong nonrel- 
ativistic shocks in a hydrogen medium, the Compton- 
scattered radiation spectrum is much harder than the 
electron bremsstrahlung 7-ray spectrum. For a consis- 
tent explanation, the synchrotron spectrum must also 
be compatible with the same electron distribution that 
makes the 7 rays. 

Differences between morphological features can help 
discriminate between leptonic and hadronic origins of 7 
rays in SNRs. Emission of 7 rays from bremsstrahlung 
and proton interactions is expected to be enhanced in 
the vicinity of dense molecular clouds, not only due 
to the denser target material, but also due to the 
increased amount of energy dissipated at the shock 
front (e.g.. lBvkov et al.ll2000l:[Aharonian fc Atovanll 19961: 
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iGabici et aO I2009D . Almost all SNRs associated 
with Fermi sources exhibit OH maser emission from 
SNR/ molecular cloud interactions (jHewitt et al.|[2009l ). 
By contrast, synchrotron X-rays and TeV 7-rays from 
Compton-scattered CMB photons require only high- 
energy electrons and a magnetic field, and in principle 
could be found in regions of small gas density. 

Recent Fermi Large Area Telescope (LAT) observa- 
tions have provided a wealth of data on SNRs. The 
First Ferm i Ca talog of Gamma Ray Sources (IFGL; 
lAbdo et al.l 12010a') lists 41 associatio ns with sources 
in Green's SNR catalog (jGreenl l2009l) . Morphological 
similarities allow definite identifications to be made in 
3 case s: IFGL J1856.1-h0122 with W44 (lAbdo et al 
20103); IFGL J1922.9-f 1411 with W51C ( Abdo et al 
203; and IFGL J0617.2-I-2233 with IC 443 ( Abdo et al 



I n the second Fermi-LAT catalog (2FGL; 



Abdo et al.ll2011al ). 3 additional SNR identifications are 
reported besides a total of 62 associations with SNRs 
and pulsar wind nebulae. More re cently, RX J 1713. 7 - 
3946 has been added to this list (|Abdo et al.ll2011bD . 
The Fermi data provide strong evidence that 7-ray emis- 
sion is made by accelerated particles in the vicinity of 
these objects. Yet it is not conclusive that the SNRs 
are accelerating the particles, as pre-existing cosmic rays 
compress ed by the outflowing r emnant could make the 
emission ([Uchivama et al.ll2010[ ). 

There have been numerous hadronic and leptonic mod- 
els produced to explain the particle acceleration and 
7-ray emission from SNRs, with varying degrees of 
complexity Ce.g..lSturner et al.lll997l: iBaring et al.llT999l: 



iZirakashvih fc AharonianI I2007t iLee et al.ll2008[ ). Here 
we focus on a simple model for electron acceleration at 
the forward SNR shock and study the evolving distribu- 
tion with the addition of radiative and adiabatic losses 
(however, note that acceleration in the rev erse shock 
could be substantial; iTelezhinskv et aLll2012[ ). The goal 
is to explain the RX J1713. 7-3946 spectrum with a purely 
leptonic model involving Compton scattering of diffuse 
target photons and a possible small contribution from 
electron bremsstrahlung. As the SNR expands into the 
GSM and decelerates due to the addition of the swept- 
up matter, the particle injection power initially increases 
in the free-expansion phase, and subsequently declines 
in the Sedov phase. The time-dependent injection effi- 
ciency is very model dependent, and ultimately rests on 
microphysical plasma processes. Here we normalize the 
power of injected nonthermal electrons to the swept-up 
power; other normalizations could employ an injection 
efficiency proportional to the rate at which particle mass 
is swept up, or an efficiency dependent on shock speed 
and compression ratio. 

In Section [2] we describe our model. In Section [31 
we use it to fit the multiwavelength SED of the SNR 
RX J1713.7-3946 (G 347.3-0.5), and show that a single- 
zone model is incapable of fitting the spectrum. A model 
with a second zone of emission regions consisting of com- 
pact knots is shown to give a good fit to the SED. We 
conclude with a summary and discussion in Section HI 

2. FORMALISM 

We make a number of common simplifying assump- 
tions. A spherically-symmetric supernova explosion finds 
itself in a homogeneous surrounding GSM with constant 



number density ncsM- The explosion is approximated 
by an expanding shell of matter that sweeps up GSM 
material. The inclusion of the swept-up mass controls 
the dynamics of the shell, and the system proceeds to 
channel directed kinetic energy into internal kinetic en- 
ergy of the shocked matter. The injection rate changes 
abruptly at the Sedov age. We purposely keep the tran- 
sition between the pre-Sedov (i.e., free expansion) and 
Sedov phases discrete to highlight interesting injection 
effects, keeping in mind that a more detailed treatment 
would have a smooth pre-Sedov and post-Sedov transi- 
tion. 

2.1. SNR Dynamics 

Here we describe a simple formalism for the dynamics 
of SNRs. A detailed h y drody namic description is given 
by iTruelove fc McKed (|199g) . The kinetic energy, E, 
of the remnant and swept up matter is conserved in an 
adiabatic blast wave, so that 



E ^ U'l^vl = i 



47r 
Mq + —mpncsMr{t) 



v^t) , (1) 



where Mq and vo are the initial remnant's mass and 
speed, respectively, rup is the proton mass, and r(t) and 
v(t) = dr/dt are respectively the radius and speed of the 
SNR as a function of time t. For simplicity, the GSM is 
assumed to be composed entirely of hydrogen. The Se- 
dov radius, r^, is defined as the radius where the mass of 
the swept-up GSM matter, 4TrmpncsMr(t)^ /3, is equal 
to the mass of the initial explosion, Mq, i.e.. 



3Afo 



1/3 



_ATrmpncsM_ 
Using Equation ^, Equation ([T]) can be rewritten 



drV 
'di) 



or 



(t) - — - «o 

^ dt ^ [l + (j./^^)3]l/2 



(2) 



(3) 



(4) 



This can be solved in the limit r <C 7*5 , giving 

r{t) = vot , 

which is known as the free expansion phase of the rem- 
nant. In the Sedov phase, r 3> Ts, giving the well-known 
behavior 

.2t, 



r{t) ^ vots 



(5) 



(6) 



Here the Sedov time is defined as 

r.. 



t. 



Vo 



The solutions ([5|) and ([6]) intersect when t 
speed of the remnant in these limits is 



v{t) 



vq 



vo {5t/2ts) 



-3/5 



r <t:rs 
r '^ Ts 



(7) 
1.84^5. The 

(8) 



The solutions for r <^ rs and r ^ rg intersect when 
t — 0.40ts- The rate at which kinetic energy is swept up 
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from the surrounding CSM is 
dE 



—rr = 2Trr^ncsMmpV^ 



(9) 



In Figure [T] we plot the radius and speed of the rem- 
nant for the exact expression, Equation Q, and the ap- 
proximate expressions, Equations (O, ([S]), and ([8]). As 
can be seen, the approximate form reproduces the exact 
behavior quite well. 

2.2. Particle Acceleration 

As the SNR shock expands into the CSM, particles will 
be accelerated at the forward shock, to which we restrict 
our treatment. We assume that the injected kinetic en- 
ergy of the nonthermal particle distribution is some frac- 
tion rj of the swept-up kinetic energy. This energy is 
swept into the shocked fluid and allows us to normal- 
ize the the nonthermal injection function, Q{'y,t) by the 
relation 



mc^ dj J Q{'j,t) = r]2TTr^ncsMmpV^ 



(10) 



where 7 is the particle's Lorentz factor and m is the par- 
ticle's mass. In order not to sweep in more energy than 
was originally available, 77 is restricted to be ^ 1, and 
the treatment is restricted to an adiabatic blast wave. 

Assuming that the injected accelerated particle distri- 
bution function can be described by a power law, then 

Qh,t)^Qoit)j-^Hir,7in2), (11) 

where H{x; a,b) — 1 ii a < x < b^ and H{x; a,b) — 



otherwise. Equation (|10p can be integrated to give 



Qo{t) = 



r'^v'^ 2TTncsMmp ij 



(g-2)(7r^-72^')-^forg^2 
X < . (12) 



[ln(72/7i)]" 



for g = 2 



Equation ^ can be inserted into Equation (IT^ and used 
to write an approximate expression for Q{"f,t), 



where 



K 



Q{-f,t)^Kf{t)j-^ 
_ v^tl 2TmcsMmp ry 



(13) 



(g-2)(7r''-72"')-' 9^2 



[ln(72/7i)] 



g = 2 



and 



fit) 
The approximation 



(r/r^ 



[1 + (r/rs) 



313/2 



J^^>~ "^ {5t/2ts)-^ t>Cts 



(14) 



(15) 



(16) 



is in accord with the asymptotes for r(t) and v{t) from 
Section 12.11 The division of the two branches of the ap- 
proximation at Cts where C w 0.74 was chosen to pro- 
duce a continuous function. The exact expression for 



/(i). Equation (J15p is compared with the approximation 
from Equation p^ in Figure [2j As can be seen, the 
approximation is quite good, with small discrepancies 
around t sa i^. 

Using the power in the swept-up CSM to normalize 
the rate of particle acceleratio n, as described above, is 
also standard for GRBs (e.g., iChiang fc Dermeilll999D 
but differs from what has been done in the past for 
SNRs. Earlier normalizations related the number of 
accelerated electrons to the numbe r of electrons swept 
into the expan ding blast wa ve (e. g., iSturner et al.lll997l : 
lRevnoldsl[T998l : baring et al.lll999f ). which leads to a con- 
siderably different time-dependence in the Sedov phase 
for the accelerated electrons {Q{-/,t) ex t^^^ instead of 

Qh,t)^t-^n 

The maximum particle energy can be calculated by 
equating the acceleration tim e with the radi ative loss 
time or the age of the remnant (|Revnoldslll998[ ). To keep 
the treatment analytic, here we assume "/max is constant 
in time. 

2.3. Particle Evolution 

When Qijjt) has been determined (Section 12.21) . the 
evolution of the particle distribution, iV(7;t) can be 
found by solving the continuity equation, 

where tesc is the escape timescale and 7 is the cool- 
ing rate. Analytic solutions to th e particle c ontinu - 
ity equation (IT71) ar e discussed in iKardashevI (I1962D: 
iBlumenthal fc Gouldl ()1970t ): and iDermer fc MenonI 
(120091, Appendix C). 

2.3.1. Solution with Radiative Losses 

For electrons, escape timescales can be long and hence 
will be neglected. The electron Lorentz factor 7 above 
which synchrotron losses dominate bremsstrahlung losses 
is 7 « 5 X 10"* (ncsM/cm-3)(10 fJ-G/B)^, but the corre- 
sponding timescale for energy loss is « 5 x 10^ yrs. This 
is much longer than the age of the remnant even for very 
dense target material, so bremsstrahlung losses can be 
safely neglected. Thus we assume electron energy losses 
are dominated by radiative losses from synchrotron and 
Thomson scattering of CMB photons, so that 



7 = 1/7 



where 



_ 4:caT[B'^/{87T) + ucMB] 
1.34 X 10-20 s-i[l + (S/3.23 ^iGf] 



(18) 
(19) 



B is the magnetic field in the remnant, and ucmb = 
4.13 X 10"^'^ erg cm^^^ is the energy density of the CMB 
at the present epoch. Klein-Nishina effects should be of 
little importance to the evolution of the electron spec- 
trum as long as the synchrotron losses dominate over 



^ Note the typographical errors in IRevnold3 IJ199SI ) and 
IBaring et all l[T999l 1 who have Q(7,t) oc t''^/^. 
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exact and approximate expressions. 




Fi g. 2. — The exact (Equation 1151 ) and approximate (Equation 
|16) ) expressions for f{t). 

Compton losses, which will be the case for B > 3 ^G. In 
this situation, the continuity equation has the solution 

^^N,{r,t) = K [ dU-ff-'^fiU) (20) 



(see Appendix X} also jKardashevI 119621 : 

IBlumenthal fc Gouldlfl970tlDermedfT998D . Here 



7i 



-f-'^ - i^it - U) ' 



(21) 



and 



i™„=max[0,t-iy-i(7-i-72-i)] . (22) 

It is instructive to look at the case q = 2, where the 
integral in Equation pop can be performed analytically. 
For t<Cts, 



7^iVe(7;0 = ^ 



(23) 



while for t > Cts and imin < Cts, 



-i^N,{r,t)^Kts 




For t > Cts and tmin > Ct 



2 Ar 1^.. ^\ _ 1 1j^ 



by] 



(25) 



Now we examine the asymptotes for this solution, 
starting with the case where t < Cts- In the limit 



t < {i^^y^ or 7 < (vt)-^, U, 



0, and 



^^N,{r.t)^K — ^t'f 



3t2 



For t > {v^y^ or 7 > [vty 

K 



^^^(^^^)-3^U 



oc t^7"\ 



(26) 



(27) 



Thus, at low 7 (7 ^ [vt) ^), the electron distribution 
will have the power-law injection index (i.e., Ne{'^]t) ex 
7^''), and at high 7 (7 ^ (i^i)~^), the electron distri- 
bution will be Nf,{j;t) oc 7"'^"''^, as seen in Equations 
((26|) and ([27|) and as expected for a cooling distribution. 
Since the electron distribution increases at different rates 
in the different regimes, however, an inflection will open 
up in the normalization of the distribution in these two 
regimes. 

This effect becomes more pronounced for t > Ctg. In 
this case, if 7 ^ (lyt)^^, 



j^N,ir,t)~Kts 



■In 



t 

Cts 



ex 7° \xi{t/Cts\2&) 



For 7 > {vtY'^ and 7 > v^'^{t - ts)~^, 



7^iV,(7;i)«^oct-i7-' 
v^t 



(29) 



Thus, at low 7 (7 <C [vt)^^), the electron distribution 
will again have the same power-law index with as the 
injection term, just as with t < Cts, and the overall nor- 
malization will increase with time, although more slowly 
than at t < Cts- At large 7, the electron distribution 
will be Ne{'j;t) ex 7"''"'"^, but the normalization will be 
decreasing with time. This will lead to an increasingly 
large gap in normalization of the cooled (7 < {iyt)~^) 
and uncooled (7 > {i^t)~^) electrons as time increases. 
This behavior can be seen in Fig. [3l Similar behavior 
can be found for q ^ 2, as demonstrated in Fig. SI where 
the integral in Equation ()20p is performed numerically. 
The parameters used in these calculations can be found 
in Tabled 

2.3.2. Solution with Adiabatic Losses 

As the remnant expands, particles lose energy due to 
adiabatic losses, since they are trapped in the expanding 
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TABLE 1 

Test model parameters. 



Parameter 


Symbol 


Model 


Blast Energy [erg] 


E 


1.0 X 10^1 


Initial Mass [Mq] 


Mo 


1.6 


Initial Veloeity [cm s~^ 


vo 


8.0 X 10« 


ICM density [cm'^] 


niCM 


1.0 


Sedov time [yr] 


ts 


303 


Magnetic field [fiG\ 


B 


10 


Cooling Constant [s^-"-] 


V 


1.4 X 10-^9 


Low energy electron cutott 


71 


10 


High energy electron cutoff 


72 


3.1 X 10» 


Electron acceleration efficiency 


'ye 


10-4 



and kad — 2/5 for r ^ r^., given r(t) and w(t) from 
Equations dS]), ([6]), and dS]). The results are only weakly 
dependent on the value of kadi ^^nd so here we will assume 
the same kad in bo th regimes , event ually taking kad = 1 
for simplicity. See iRevnolda (12008^) for a discussion on 
the dependence of the shock radius with time for core 
collapse and type la progenitors, which has implications 
for kad- 

If adiabatic cooling dominates and radiative cooling is 
negligible, then the solution to the continuity Equation 
gives 



N{r,t) = KtJ^ ^ 



^i^kaa 



liU) (31) 




Fig. 3. — Evolution of electron distribution with radiative losses 
only for q = 2. Parameters are given in Tabic [T] 
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Fig. 4. — Same as Figure[3]for q = 2.5. Parameters are given in 
Table [1] 



SNR. The loss rate from this process is 
-1 dV I kadi 



'^ i'dt V 



t 



(30) 



(e.g., iGouldl [r975[ ) where the remnant's volume \s V (x 
r(tY ■ Details of the expansion control the adiabatic co- 
efficient kad, and in principle one can imagine a shell 
that contracts in width while expanding outward to give 
no volum e change, so kad ~^ for this peculiar system 
(see, e.g. iTruelove fc McKed [19991 for a hydrodynamic 
description). In our approximations, kad = 1 for r <C r^. 



where 



l/fead 



(32) 



The integral can be performed analytically. If t < Cts , 



N{r,t) 



Kt^j-'J 



tl{qkad + 2) 
If t > Cts and 7t/72 < ts, 



qkad+2 



(33) 



t' 



'Ct. 



tl{qkad + 2) W t 



qkad+2 



N{-/;t) = K-f-'i 
(qkc,d+'2)/Kd\ 



2t, 



?>{qkad - 1) 

If t > ts and 7^/72 > ts, 
2K-f-n 




N{r,t) 



5{qkad - 1) 



1-1^ 

72 



(qkad-i)/kad 



.(35) 



2.3.3. Solution with Radiative and Adiabatic Losses 

If adiabatic and cooling losses are important, then the 
cooling rate will have the form 



7 = 1/7 + kad 



7 



(36) 



In this case, if kad — 1, the continuity Equation has the 
solution (see Appendix IBJ) 



K 



^2-9 



YN{r,t) = - dt.t.-fr" fit,) (37) 



where 



7^ = [u{{ity^-'^Ht/u))Y 



W ( —e^/(-'r''^' I 72 i^ 



(38) 



(39) 



\72J^t 

and W{x) is the Lambert W function, defined by 

X = I^(x)e^(") (40) 
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(e.g.. iCorless et aLlll996[ ). The integral can be done an- 
alytically if g = 2. In this case, for t < tg, 



I'Nir.t) 



K 

ITTf 

For t > Cts and i„i„ < Ct^, 



y ''minj 



j^N{r,t)^j 



1 Of 



For t > Cts and U 



>Cts 



i'N{r,t) 



2Kts 
5t 



yi t^min) 



(41) 



'42) 



(43) 



The asymptotes of this solution can also be found. We 
begin with the case oi t < Ctg. For t ^ {t^l)~^ or 
7 <C {vt)~^, the argument in the Lambert function in 
Equation (1391) goes to infinity, so the Lambert function 
goes to infinity, and tmin — >■ 0. Then, 



7^iV(7;i)«^« t'l'- 



(44) 



For t ^ (1^7) ^ or 7 ^ (i/t) ^, a Taylor expansion of 
the exponential and Lambert function in Equation (j39p 
gives tmin -^t- (j^7)~^ Then 



7^N{r,t) 



Kt^ 



IGt^v^ 



ex t^-^~^ 



(45) 



These results are quite similar to the asymptotes for the 
radiative-only case (Section 12.3. ip . The results for t > 
Cts are, however, somewhat different. In this case, for 
i< (j^7)~^ or 7< (lyty^, 



^^N{r,t)^Kts 



3U 
20t 



2Kt, 



oc<%o 



For t > {v-iY^ or 7 > {vty 



51^17 



(46) 



(47) 



Here, the behavior for large 7 is essentially the same as 
the radiative-only case. However, the addition of adia- 
batic losses means that the electrons at low 7 will ap- 
proach a constant value, rather than increasing logarith- 
mically without bound. This will remove the inflection 
seen in the radiative-only case. An example of this can 
be seen in Figure [SJ where it is clear the inflection in 
the radiative losses-only solution (Section 12.3. ip is not 
seen. Results are very similar for other values of kad, 
particularly kad = 0.4, as one would expect in the Sedov 
phase. 

2.4. Spectral Energy Distribution 

Once the electron distribution has been determined, 
as above, the spectral energy distribution (SED) from 
the SNR can be calculated. For electrons in a randomly- 
oriented magnetic field B, the lyF^ synchrotron fiux, f^^^ 
is given by 



f: 



AirhdP 



dj N,ir,t) Rix) , (48) 




Fig. 5. — Evolution of electron distribution for q=2 including 
adiabatic and radiative cooling. Parameters are given in Table [T] 

where e = 4.8 x 10"^" esu is the elementary charge, e = 
hv l(meC^) is the dimensionless observed photon energy, 
h = 6.63 X 10~^^ erg-s is Planck's constant, d is the 
distance to the SNR, 



47rmgC^ 
3eBhj^ 



R{x) = 



de sin 61 



x/ sin 9 



(49) 



dt K^,^{t) (50) 



(jCrusius fc Schlickeiseii I1986D . and K^/'j,{t) is the 
modified Bessel function of or der 5/3. Approximate ex- 
pressions for R{x) are given bv'Zirakashvili fc AharonianI 
(,2007.): Finkc, Dormer, feBottchcr (200^^ and 
fjoshi fc Bottcheri (|2011D . 

The electrons will also Compton-scatter external pho- 
ton sources, such as the CMB or other intergalactic 
sources. The flux from external Compton (EC) scatter- 
ing a blackbody photon source with total energy density 
utot and dimensionless temperature Q — ksT / (raeC?) is 
given by 



/f" = 



3TOeCCTTe^ 15Mfo( 



(51) 



de* 



exp(e,/9) — 1 

d7 7-2iVe(7;i) i^c(e,7,e*) 



where ctt — 6.65 x 10^^^ cm^ is the Thomson cross 
section. For the CMB at the present epoch, 8 = 
4.58 X 10"^° and utot = 4.2 x 10"^^ erg cm"^ are the 
dimensionless CMB temperature and energy density at 
the present epoch. The integration limits are given by 



_ 1 



1 



e*e 



1 



1 
e*e 



(52) 



Iniax = ff (e* - e) , (53) 

e* — e 



and the function 



2,ln,+ (l+2,)(l-,) + iJi;Hf^(l-,) 
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where 



X H 



e/7 



1 



9;77^>i 



and r. 



4e,7 



(54) 



(55) 



re(l - e/7) 

ponesi 119681: IBlumenthal fc GouldlllQTOl) . The electrons 
can also Compton-scatter the synchrotron photons pro- 
duced by the same electron population (known as syn- 
chrotron self-Compton or SSC), which is given by 



/, 



SSC 



9 crxe^ 
16 Trr^ 



de. 



J^syn 



d^^Y^Fc{q,r) (56) 

(e.g-. lFinke et al.ll2Q08[ ). This mechanism is usually neg- 
ligible for large, diffuse remnants, but will play a roll in 
Section [3Jl 

The nonthermal electrons from the remnant will inter- 
act with the cold ions (assumed to be protons) in the 
surrounding CSM to make bremsstrahlung (or free-free 
radiation) with flux given by 



/, 



47rd2 



d7iVe(7;i)^^^(^;7), 



(57) 



where the bremsstrahlung cross section is written as 

dcTff^eZ 



de 



-(e;7) 



AZ'^afrl 



1 



2y 



In 



2-fhj 



(58) 



(jBlumenthal fc Gouldlll970l ). Z is the effective charge of 
the cold ions, a/ « (137)"-^ is the fine structure constant, 
r-g = 2.82 X 10"^'^ cm is the classical electron radius, and 
y = l- e/7. 

The synchrotron and Compton emission can be seen 
in Figs. inland [71 corresponding to the respective electron 
distributions seen in Figs. [31 and |5l These SEDs assume 
emission from an SNR at d = 1 kpc. In the Sedov phase, 
the feature resulting from the different evolution above 
and below the cooling break is clearly seen for the case 
when adiabatic losses are neglected (Figure (5]) , but this 
feature is not as pronounced when adiabatic losses are 
taken into account (Figure [7]). 

3. APPLICATION TO THE REMNANT RX J1713. 7-3946 

We apply our results to SNR RX J1713. 7-3946 
(G 347.3—0.5), which is thought to be the remnant of 
a "g uest star" observed by Chinese astronomers in 393 
CE (jWanget all [1991 . This fixes it age at ^ 1620 
yr, so SNR RX J1713. 7-3 946 is likely to be well into 
the Sedov phase, although i Fukui et al.l (|2003D argue in- 
stead that it is still in the free expansion stage, for 
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CSM densities ng s M < 0.01 cm^^ at sa 1 kpc dis- 
tance. iSlane et al.l ( 19991 ) associated the source with a 
nearby molecular cloud, giving a distance to the source of 
^ 6 kpc. However, n ewer CO observations found molec- 
ular gas at ~ 1 kpc, (jFukui et al.ll2003l : iMoriguchi et al.l 
i2003). Absorbing column densities from X-ray ob- 
servations strengthen the ^ 1 kpc d istance estimate 
(jKoo et al.ll2004 iCassam-Chenai et al.ll2004D . making it 
the most likely one. The detection of an X-ray point 
source, thought to be a left-over neutron star, implies 
the remnant is the r esult of a core-collapse supernova 
(jLazendic et al.ll2003D . 

The X-ray spectrum of RX J1713. 7-3946 appears 
completely domina ted by nonthermal emission (e.g., 
iTanaka et al.l [20081 ) . The lack of thermal X-ray lines is 
taken as evidence that npsM 5, 0.2 cm~'^ (|Slane et al.l 
[T99I lEllison et al.llMot lAbdo et al.ll2011bD . Without a 
dense target for cosmic ray protons, 7r° decay is prob- 
ably not a significant contributor to th e 7-ray spec- 
trum of this remnant (l EUison et al.|[20lol ). and electron 
bremsstrahlung is weak in comparison with the Comp- 
ton 7-ray emission. But note the different energy ranges 
of electrons that radiate into the 7-ray band: electrons 
with Lorentz factor 7 -- 10*5 (-- 10^) scatter CMB pho- 
tons to GeV (10 TeV) ranges, and electrons with 7 ~ 10^ 
(~ 10^) make GeV (TeV) bremsstrahlung. 
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HESS observations show that the X-ray image 
of RX J1713. 7-3946 a nd the VHE 7-rays a r e spa - 
tially well-corre l ated (jAharonian et al.l 12004 I2006D . 
lUchivama at alJ ()2007l ) reported X-ray variability on a 
year timescale in a few small (arcsecond scale) hotspots 
of RX J1713. 7-3946. If this reflects radiative variabil- 
ity of the nonthermal electrons, large magnetic fields are 
needed (B ^ 1 mG), and the implied small number of 
electrons in such a strong magnetic field could not pro- 
duce the TeV emission. But the appearance of thin radio- 
emitting rims in some SNRs could mean that the emis- 
sion zone is compact and variability is due to shell expan- 
sion or compression as it encounters d ilute or dense GSM 
(|Revnoldsll2010HRevnolds et al.ll201lD . Furthermore, the 
existence of knots which do not seem to be variable on 
such short timescales indicates knots may exist with sig- 
nificantly lower magnetic field s. Ba sed on the HESS 
observations, Aharonian et"al] (j2006() concluded that a 
leptonic model was unlikely to fit the broadband SED, 
and a hadronic origin was favored for the 7 rays from 
RX J1713. 7-3946. The variable X-ray filaments cannot , 
however, explain global TeV emission (jButt et al.ll2008D . 

With the arriva l of the first epoch Fermi data from 
RX J1713. 7-3946 (lAbdo et al.l [2OIIB . different models 
can be tested better. This seems to make it a good 
time to revisit leptonic models, which are natural for 
RX J1713. 7-3946 given the close spatial correlation be- 
tween the X-ray and TeV 7 rays. 

3.1. Single Zone Model Fit 

The integrate d broadband SEP o f RX J1713. 7-3946 is 
shown in Fig. [5] iPorter et al.l (|2006f ) have shown that the 
interstellar infrared radiation field (IIRF) can be a signifi- 
cant photon source for Gompton scattering. It is strongly 
dependent on the position in t he Galaxy, and close to 
the Galactic center their model (jMoskalenko et aLll2006l) 
gives the IIRF energy density greater than that of the 
GMB. At a Galactic longitude of £ = 347.3° and a 
Galactic latitude of 6 = —0.5°, this remnant is nearly 
along the line of sight of the Galactic center. This makes 
the intensity of the IIRF strongly dependent on its dis- 
tance from the Earth (and therefore the Galactic cen- 
ter). The consensus, based on associations of molecu- 
lar clouds and absorption of X-rays (as discussed above, 
in Section [Sj seems to be that the RX J1713.7-3946is 
at d = 1 kpc from Earth. Therefore in our modeling 
we use this distance, and an IIRF intensity consistent 
with this distance (or 7.5 kp c from the Galactic cen- 
ter) from IPorter et al.l ()2006D . modeled as a blackbody 
with temperature T = 30 K and total energy density 
u = 4.8 X 10" ^^ erg cm'^ = 0.30 eV cm-^. Note that 
ILi et al.l (J2011D modeled the sourc e using the IIRF a t a 
level near what one would expect (" Porter et al.ll200"6l ) if 
RX J1713.7-3946 was d = 6 kpc from the Earth. 

The models in Figure [5] include adiabatic and radia- 
tive losses from both the GMB and IIRF, and apply to 
the integrated emission over the entire remnant. In our 
modeling we have added this radiation field as an addi- 
tional term in equation (1191) to take it into account in 
the evolution of the SNR. We take the age of the rem- 
nant to be t = 1620 yr. A best fit to the SED, given 
reasonable parameter constraints on age, ICM density, 
initial mass, and blast energy, is shown in the figure as 
the black curve. This model includes emission from syn- 
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Fig. 8.— The integrated broadband S EP of RX J1713 . 7-394 6 
with data from ATCA (green triangles; I Aharonian et al.' 1200(3), 
Suzaku (blue circle s; [^naka ct al. 2008), LAT (black diamonds; 
lAbdo et al.l l2 011bD , and HESS (red diamonds; .Aharonian et all 
2006, 2007). The solid curves show the model fits for different 
magnetic fields, as labeled. The dashed and dotted curves show 
the Compton-scattered CMB and IIRF components, respectively, 
for Model 1. 



chrotron, Compton-scattering of GMB and IIRF pho- 
tons, and bremsstrahlung, which can be seen as a small 
bump at ~ 10^^ Hz. It provides a good fit to the radio 
and X-ray data, but does not adequately reproduce the 
LAT and lower-energy HESS measurements. 

The model parameters are shown in Table [21 The pa- 
rameters which are constrained by the SED fit are B, q, 
72, and rj^. We also display in this Figure models with 
5 times larger and small B. As the magnetic field in- 
creases, the overall synchrotron flux increases, and the 
cooling electron Lorentz factor {vt)~^ decreases. The 
magnetic field will only affect the Compton-scattered 
emission through its effects on the electron distribution; 
thus the overall Gompton-scattered fiux will not increase, 
and indeed will decrease above the the cooling break, 
which is lower for higher B. 

In Fig. [9] we explore variations in uq , which also cor- 
respond to variations in Mq, assuming E is held con- 
stant, as given in Equation ([T]). The Sedov radius 

Ts (X Vq as given in Equation ([2|) , and thus i^ oc Uq 

(Equation [7]) and K o^ Vq (Equation [14]). Keeping 
this in mind, i > i^ so that below the cooling break, 

Ne{T,t) oc 1/3 + (5/2) ln(wo/^) from Equation ^ and 
above the break, Ne{j;t) ex Vq (Equation [29]). This 
is reflected in the SED, as seen in Fig. |9l where below 
the break the emission increases with uq, while above the 
break, the emission is independent of vq. Also note that 
it is unlikely for the initial ejecta mass to be as low as 
found in Model 4, demonstrating the limited usefulness 
in varying vq to obtain a good fit to RX J1713. 7-3946. 

3.2. Multi-Zone model 

As discussed above in Section [S] the discovery of vari- 
ab le X-ray filament s (or knots) within the SNR structure 
bv lUchivama et ahl (2007) indicates that the single-zone 
fit is inadequate to explain the overall SED. However, 
the filaments themselves could contribute a significant 
amount to the 7-ray emission from the source. 

Smaller knots emitting synchrotron, SSG, and Comp- 
tonized CMB and IIRF radiation were added to Model 
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TABLE 2 
RX J1713. 7-3946 Model Parameters 



Parameter 



Symbol Model 1 (baseline) Model 2 



Model 3 



Model 4 



Model 5 



Blast Energy [erg] 
Initial Mass [Mq] 
Initial Velocity [cm s" 
ICM density [cm"^] 
Sedov time [yr] 



E 

Mo 

vo 

niCM 

ts 



1.6 X lO^i 

1.6 
1.0 X 10^ 

0.2 

420 



1.6 X 10^1 

1.6 
1.0 X 10^ 

0.2 

420 



1.6 X 10">i 

1.6 
1.0 X 10^ 

0.2 

420 



1.6 X 10=1 

6.4 

5 X 10** 

0.2 

1300 



1.6 X IQSi 

0.4 
2.0 X 10^ 

0.2 

132 



Magnetic field [/jG] 


B 


12 


60 


2.4 


12 


12 


Cooling Constant [s"^] 


V 


2.2 X 10"^'' 


4.7 X 10^1* 


3.7 X 10-2" 


2.2 X lO-i'* 


2.2 X 10-19 


Cooling electron Lorentz factor 


(yt)-^ 


9.1 X lO'^ 


4.2 X lO'' 


5.4 X 10* 


9.1 X 10^ 


9.1 X lO'^ 



Low energy electron cutoff 71 10 10 10 

High energy electron cutoff 72 3.1x10* 3.1x10* 3.1x10* 

Injection spectral index q 2.1 2.1 2.1 

Electron acceleration efficiency r^e 5.0 X 10"^ 5.0 x 10^^ 5.0 X 10~^ 



10 

3.1 X 10* 

2.1 
5.0 X 10-5 



10 
3.1 X 10* 

2.1 
5.0 X 10-5 



1^ lii^ 111^ iiii^ 11^ lii^ 111^ iiii^ 11^ iii^ 111^ iiii^ 11^ lii^ 111^ iiii^ 11^ 111^ iiii^ iiii^ 11^ 111^ 
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Fig. 10.— Multi-zone model fit to RX J1713. 7-3946. Curves 
show the total emission from the knots and overall shock combined 
(black solid curve) as well as synchrotron emission (dot-dashed 
curves), Compton-scattered CMB (dashed curves), Compton- 
scattered IIRF (dotted curves), and SSC (double dot-dashed 
curves) from the overall shock and knots. 



1 , as seen in Figure [TOl The much smaller volume of 
these knots results in large synchrotron energy densi- 
ties in the knots, with strong SSC emission at GeV en- 
ergies. This fit has the number of zones taken to be 
^knots = 100, with each zone having Bknota = 16 A*G, 
radii Rknot = 1 mpc, and an electron distribution that 
spans from ^knot,i = 10 to jknot,2 = 1-4 x 10* with a 



break at ^knot,brk = 4.7 x lO"^ with Ne,knot{l) oc 7 ^-^ for 

7 < lknot,brk and Ne,knoth) OC 7-^'^ for 7 > Jknot.brk- 

As can be seen in Figure ^0\ this reproduces the SED 
well, and makes interesting predictions. 

The synchrotron component is dominated by the large 
first zone that effectively represents the entire remnant, 
which also makes the bulk of the TeV radiation. Emis- 
sion > 1 TeV is dominated by the Compton-scattered 
CMB of the remnant as a whole, while in the range in 
the joint LAT/HESS window from < 1 TeV the 7 rays 
arise from the SSC component in the knots. The angu- 
lar resolution of the LAT is generally worse than 0.1°. 
At a distance of 1 kpc, the 1 mpc knots will have an 
angular radius of 0.2" and thus cannot be resolved with 
the LAT. CTA will have an angular resolution of ^ 1' 
(jCTA Consortium 2010) and will not be able to distin- 
guish the variab l e and non- variable X-ray knots seen by 
lUchivama et al.l (|2007[ ) either, even if they radiate in 7 
rays. However, if the low and high energy 7-rays come 
from different components, maps of RX J1713. 7-3946 
made with CTA may be different at lower (< 1 TeV) 
and higher (> 1 TeV) energies, with the higher energy 
maps being more in agreement with X-ray ones. This 
may allow this multi-zone model to be tested. 

The knots contribute ~ 10% to the X-ray emis- 
sion of the remnant, c onsistent with observations from 
lUchivama et al.l ()200 3D. They are also much lower than 
the va lues inferred from variability by lUchivama et al.l 
(|2007f l. However, there seem to be many knots which 
are not variable, which could reflect a lower magnetic 
field. 

4. DISCUSSION AND CONCLUSIONS 

The SNR RX J1713. 7-3946 occupies an important 
place in 7-ray studies of supernova remnants. TeV emis- 
sion from the SNR RX J1713.7-3946 wa s first detected 
with the CANGAROO experiment (jMuraishi et al.l 
2OOOI). Based on further CANGAROO observations. 



Enomoto et al.l ()2002 f) claimed that a standard leptonic 



synchrotron/EC-CMB model did not fit these data, in- 
cluding the EGRET upper hmit. iReimer fc PohJI ()2002[ ) 
argued that EGRET upper limits rule out a hadronic ori- 
gin, but diffusion of high-energy part icles upstream of the 
shock c an harden nuclear emission (jMalkov fc Diamondl 
I2006f) . lAharonian et al.l (|2004[ ) produced the first re- 
solved 7-ray image of an SNR by HESS. Further HESS 
observations found that the X-ray and VHE 7-rays 
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were spatial ly well -correlated (jAharonian et alj I2006D . 
iPorter et alj (|2006D found, however, that Compton- 
scattered Galactic background photons, in addition to 
CMB photons, could help to explain the RX J1713.7- 
3946 VHE emission in leptonic models. Still fur- 
ther HESS observations detected the remnant o ut to 
~ 100 TeV (jAharonian et al.ll2007[) . iLi et all (|20TT1 ) pro- 
vide a good fit to the full SED including the LAT spec- 
trum with a model similar to IPorter et al.l (|2006[ ). in- 
cluding Compton-scattering of interstellar infrared pho- 
tons. As discussed above in Section 13. 1) they assumed 
the source was at a distance of 6 kpc from us, closer to 
the Galactic center where the IIRF is much more intense. 
However, we think the molecular cloud and X-ray ab- 
sorption evidence points to RX J1713. 7-3946 most likely 
being at d = 1 kpc. This emphasizes the crucial im- 
portance of an accurate distance measurement to SNR 
modeling. 

New data from Fermi (jAbdo et al.l l2011bD , in ad- 
dition to multiwavelength measurements at radio, X- 
ray, and TeV energies, reveal the bolometric SED 
of SNR RX J1713.7-3946 with unprecedented de- 
tail. The joint Fermi-LAT/H ESS data favor models 
(IPorter et al. 2006; Berezhko & Volk 200l; lEUison et all 
l2010t iZirakashvili fc Aharonian.2010. ) where the 7 rays 
have a leptonic rather than a hadronic origin. This 
conclusion follows rather forcefully if the injection spec- 
tral index of the particles — protons or electrons — is 
softer than g = 2, as expected in the linear first-order 
Fermi acceleration theory. Nonlinear effects m ay mod- 
ify the injection index (e.g., iBlasi et al]|2005[) . Indeed, 
lYamazaki et al.l ()20Q9f ) point out that nonlinear effects 
could harden the emission in the LAT energy range, mak- 
ing it nearly impossib le to d istinguish between leptonic 
and hadronic origins. Ilnoue et al. (2011) find that hard 
7-ray spectra can be generated from 7r° decay if the GSM 
has inhomogeneities and is "clumpy" . 

The se results also test the conclusions of iFukui et al.l 
(|2011| ) based on a comparison of the TeV and GO and 
H I morphology. They interpreted the good correlation 
between the two bands as being a strong signature for a 
hadronic origin of the 7 rays, since cosmic ray protons 
would react with the molecular cloud hadrons. That in- 
terpretation is not unique to hadrons, however, as shocks 
in molecular clouds would enhance electron acceleration 
and leptonic 7-ray emission. Observations with HAWG 
and GTA, and longer exposures with the Fermi-Ij AT , will 
show the fraction of radiation from clumps at different 
energies, and will help clarify the issue. Here we con- 
sider broadband spectral modeling by following electron 
injection and evolution. 

The complicated GSM distribution in any realistic 
SNR environment is quite different from the assumption 
of a homogeneous medium, but within this approxima- 
tion, we reconsidered particle injection, and found that 
the assumption that the injection power is proportional 
to the rate at which kinetic energy is swept downstream 
of an adiabatic blast wave yields interesting structure in 
the particle injection distribution that is cooling. The 
addition of adiabatic losses significantly smooths these 
effects, but in either case, fitting of the RX J 1713. 7-3946 
data with a single-zone synchrotron/Gompton-scattered 
model did not give a perfect fit. 



The addition of kn ots, as in the two-zone models of 
lAtovan et al.l (|2000al lH) applied to Gas A, introduce in- 
teresting effects on electrons escaping downstream into a 
region of different magnetic field. In a two-zone model, 
particles may be accelerated in smaller kn ots, and diffu se 
into a larger zone. The model of lAtovan et al.l (|2000ai rbl) 
is justified by the small knots seen in the 6.3 cm Very 
Large Array (VLA) image of the remnant Gas A. The 
VLA or another high angular resolution radio telescope 
has not yet observed RX J 1713. 7-3946, although we jus- 
tify this complica tion from knots obse rv ed from this 
source in X-rays (jUchivama et al.l[2007| l. iBvkov et al.l 
(|2008D showed that the structure and variability in these 
X-ray images can be reproduced with a steady electron 
population in a random magnetic field. Similar variable 
(on ~ 4 y ear timescales) X-ray knots h ave been found 
in Gas A (jUchiyama fc AharonianI I2OO8I ) . also implying 
large fields (S ~ 1 mG) similar to RX J1713.7-3946, as- 
suming the variability is due to rad iative cooling. Gas A 
is the result of a lib supernova (jKrause et al.l I2008D , 
and RX J1713. 7-39 46 is probably the re sult of a core- 
collapse supernova (jLazendic et al.l [20030 . so they seem 
to be of similar type, although Gas A is much younger 



than RX J1713. 7-3946, ^ 300 years (Fe sen et a 
versus ^ 1600 years. Note that .Yamazaki et al 



2006) 



Wm 



have also applied two zone leptonic and hadronic models 
to RX J1713. 7-3946, finding that one-zone models could 
not explain the VHE 7-ray spectrum and LAT upper 
li mits available at t h e time . 

iKatz &: WaxmanI ()2008n use radio observations of 
SNRs in nearby galaxies to put a lower limit on the ra- 
tio of accelerated electrons to protons. Assuming this 
ratio is approximately the same for all SNRs, they find 
that a hadronic explanation fo r the 7-ray emission from 
RX J1713. 7-3946 is unlikely. lYuan et al . (2011a) have 
fit the broadband SED of RX J1713. 7-3946 with three 
models, consisting of leptonic, hadronic, and hybrid lep- 
tonic/hadronic emission. They found their hadronic 
model provided the best fit, but it also had the great- 
est number of poorly-constrained free parameters. Be- 
cause of this, along with requiring an unrealistically large 
amount of energy put into nonthermal protons, they 
concluded that they could not distinguish between their 
three scenarios. 

Based on spectral modeling of the broadband SED of 
the Tycho SNR, particularly the shape of the Fermi- 
LAT and VERITAS spectra, it has been suggested 
that only hadronic emission, and not leptonic emis- 
sion, can be the source of 7 rays from this object 
(|Morlino fc Gapriolil 12011 IGiordano et~alll20lll) . How- 
ever, this has been called into question by a two-zone 
model (JAtoyan & Dermer 2011). Final conclusions re- 
garding cosmic-ray proton/ion acceleration in Tycho rest 
on the spectral shape below sa 400 MeV. A two-zone lep- 
tonic model for RX J1713. 7-3946, as we have seen here, 
avoids any need for cosmic-ray proton acceleration. Note 
that this does not preclude cosmic-ray proton accelera- 
tion either, simply that protons do not contribute signif- 
icantly to the emitted electromagnetic radiation. 

In conclusion, we have described a simple model for 
the time evolution of SNR emission. This includes an as- 
sumption that particle acceleration efficiency is propor- 
tional to the power swept into the expanding blast wave. 
Effects of radiative and adiabatic cooling on the evolv- 
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ing particle distribution, and emission from synchrotron, 
bremsstrahlung, and Compton-scattering processes were 
taken into account. In doing this, we have made a num- 
ber of simplifying assumptions. We have assumed the 
CSM density is constant, which may be less likely for 
remnants of core-collapse rather than Type la super- 
novae. We have assumed the magnetic field strength 
and power-law injection index do not vary with time, 
and have used a simple Sedov solution neglecting re- 
verse shocks. Hadronic emission processes were neglected 
in this study. We have applied our evolution model to 
RX J1713. 7-3946, and showed that a single-zone model 
cannot reproduce its SED if it is at a distance d = I kpc. 
The addition of a second zone consisting of compact 



knots gives an acceptable fit and makes interesting ra- 
dio and 7-ray predictions that should be testable in the 
near future. 
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APPENDIX 
A. SOLUTION TO THE CONTINUITY EQUATION FOR RADIATIVE COOLING ONLY 
We wish to solve the continuity equation. Equation (J17p . for tesc -^ 00 and 7 = —v^'^: 

^--^b'Nir,t)]-Qh,t). 



This has the solution 



N{r,t) 



t /'72 

dti / dji G{ji,ti,j;t) Q{'yi,ti) 

•'71 



(Al) 



(A2) 



where 71 and 72 are the respective lower and upper limits on the injected electrons' Lorentz factors, G'(7i,ii,7; t) is 
the Green's function which satisfies the equation 



a^--^[7'G]=<5(7-70'5(t-t.) 



(A3) 



and d{x) is the standard Dirac delta function. A single particle injected with Lorentz factor 7^ at time t^ will, at some 
later time t, have a Lorentz factor ^(t) which satisfies the equation 



dl 2 



Equation (|A4p can be solved, 



(A4) 



7(t) 



7r' + Ht - U) 



and thus the solution to Equation (|A3p is 
which can be rewritten as 



(A5) 



(A6) 



7 
Inserting this into Equation (jA2l) . one can perform the integral over 7^ to get 

7 Jt„i„ -■ L*- 

where 

imm = niax[0,i- i^"\7"^ -72"^)] . 
The lower limit tmm comes about because particles are injected only with 7^ < 72. 



(A7) 



(A8) 



(A9) 
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B. SOLUTION TO THE CONTINUITY EQUATION FOR RADIATIVE AND ADIABATIC COOLING 
We now wish to solve the continuity equation, Equation ([TTl) for radiative and adiabatic cooling, i.e., 

We can follow the same procedure as in Appendix [Xj Equation (JB1[) can be solved for 7(t) to give 



7(t) = |ife-[(74-)-i + z.T(t,t,)]} ' 



(|Gupta et al.ll2006[) where 



T{t,U 



Ht/u) 



ti"^')/{l~kad) kad^l 



In this case Equation (jB2|) implies the Green's function which satisfies Equation (|A3P is 



G{-f„u,r,t) 



1 



j^t^^..t^^^[{^t)-^-,yT{t,U)Y 



1 



It 



i^^-[(^iA..)-i_^T(i,tO] 



(Bl) 

(B2) 
(B3) 

(B4) 



Substituting this into Equation (jA2p and performing the integral over 7^ with the help of the Dirac (5-function gives 



N{r,t) 



K 



tkad 



dU t^-" f{U) 



The lower limit tmm can be found from the constraint that 

1 



U[{^t)-^-vT{t,t,)] 



2-q 



For kad = 1, when solved for U this constraint gives 



< 72 



^min, — 



\l2vt 



/hvA 



n -1 






72 V 



(B5) 



(B6) 



(B7) 



where W{x) is the Lambert W function (jCorless et al.l [19961 ) . For general values of kad. Equation (jB6p does not have 
a simple analytic solution, and it is solved numerically for tmin- 

If kad = 0, i.e., there are no adiabatic losses, then Equation (IB2p will reduce to Equation (jASp . leadin g to the 
radiative losses-only solution. On the other hand, if radiative losses are negligible, i.e. i' — ?► 0, then Equation (jB2p will 
reduce to 7(i) = 7i(ii/i)'^'"* which leads to the solution in Section r2. 3. 2 1 
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